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We study the B 4>Ks and B ri'Ks decays in MSSM by calculating hadronic matrix elements 
of operators with QCD factorization approach and including neutral Higgs boson (NHB) contribu- 
tions. We calculate the Wilson coefficients of operators including the new operators which are 
induced by NHB penguins at LO using the MIA with double insertions. We calculate the as order 
hadronic matrix elements of the new operators for B 4>Ks and B r]' Ks- It is shown that the 
recent experimental results on the time-dependent CP asymmetries in _B ~+ (j>Ks and B — > r]' Ks, 
S^K is negative and Sr,'K is positive, which can not be explained in SM, can be explained in MSSM if 
there are flavor non-diagonal squark mass matrix elements of 2nd and 3rd generations whose size sat- 
isfies all relevant constraints from known experiments {B — > Xs"f,Bs jj.'^fj,~,B XsfJ,^ j-L~ , B 
Xsg, AMs, etc.). In particular, we find that one can explain the experimental results with a flavor 
non-diagonal mass insertion of chirality LL or LR or RR when Os corrections of hadronic matrix 
elements of operators are included, in contrast with the claim in the literature. At the same time, 
the branching ratios for the two decays can also be in agreement with experimental measurements. 



I. INTRODUCTION 



The measurements of the time dependent CP asymmetry S'j/^r- in B ^ J/ijjKs have established the presence of 
CP violation in neutral B meson decays and the measured valueu] 

world — avc 

= 0.734 ± 0.054. (1) 

is in agreement with the prediction in the standard model (SM). Recently, various measurements of CP violation in 
B factory experiments have attracted much interest. Among them 

S^Ks = -0.39 ± 0.41, 2002 World - average 

S^i^^ = -0.15 ±0.33, 2003 World - average (2) 

is especially interesting since it deviates greatly from the SM expectation 

S^Ks = sin(2/3(0Xs)) = sm{2(3iJ/^PKs)) + 0{X'^ (3) 

where A ~ 0.2 appears in Wolfenstein's parameterization of the CKM matrix. Obviously, the impact of these ex- 
perimental results on the validity of CKM and SM is currently statistically limited. However, they have attracted 
much interest in searching for new physics 0j S S S 13 it has been shown that the deviation can be understood 
without contradicting the smallness of the SUSY effect on _B ^ J/^pKs in the minimal supersymmetric standard 
model (MSSM) 0,|3. 

Another experimental result, which is worth to notice, is of the time dependent CP asymmetry Srj'Ks in -B — > 

Srj'Ks = 0.02 ±0.34 ±0.03 BaBar 

= 0.43 ±0.27 ±0.05 Belle (4) 

which deviates sizably from the SM expectation. Although both the asymmetries Scj,K and S^/'k are smaller than 
the SM value, S^k is negative and Sri'K is positive. Because the quark subprocess b — > sss contributes to both 
B — > (j}Ks and B — > r]' Ks decays one should simultaneously explain the experimental data in a model with the same 
parameters. It has been done in Ref.^^l ^ model- independent way in the supersymmetric (SUSY) framework. In 
Ref.|lOl| the analysis is carried out using the naive factorization to calculate hadronic matrix elements of operators 
and the neutral Higgs boson (NHB) contributions are not included. As we have shown in a letter 8] that both the 
branching ratio (Br) and CP asymmetry are significantly dependent of the corrections of hadronic matrix elements 
and NHB contributions are important in MSSM with middle and large tan (3 (say, > 8). In the paper we shall perform 
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a detailed analysis of SiI>Ks Srj'K as well as Br in MSSM including NHB contributions and the as corrections of 
hadronic matrix elements. 

We need to have new CP violation sources in addition to that of CKM matrix in order to explain the deviations 
of S^Ks and Sri'K from SM. There are new sources of flavor and CP violation in MSSM. Besides the CKM matrix, 
the 6x6 squark mass matrices are generally not diagonal in flavor (generation) indices in the super-CKM basis in 
which superfields are rotated in such a way that the mass matrices of the quark field components of the superfields are 
diagonal. This rotation non-alignment in the quark and squark sectors can induce large flavor off-diagonal couplings 
such as the coupling of gluino to the quark and squark which belong to different generations. These couplings can 
be complex and consequently can induce CP violation in flavor changing neutral currents (FCNC). It is well-known 
that the effects of the primed counterparts of usual operators are suppressed by rris/mi, and consequently negligible in 
SM because they have the opposite chirality. However, in MSSM their effects can be significant, since the flavor non- 
diagonal squark mass matrix elements are free parameters which are only subjective to constraints from experiments. 

For the 6 — > s transition, besides the SM contribution, there are mainly two new contributions arising from the 
QCD and chromomagnetic penguins and neutral Higgs boson (NHB) penguins with the gluino and squark propagated 
in the loop in MSSM* . The relevant Wilson coefficients at the mw scale have been calculated by using vertex mixing 
method in Ref . J^] . There is the another method, " mass insertion approximation" (MIA) , which works in flavor 
diagonal gaugino couplings gqq and diagonal quark mass matrices with all the flavor changes rested on the off-diagonal 
sfermion propagators. The MIA can be obtained in VM through Taylor expansion of nearly degenerate squark masses 
rriq- around the common squark mass rrig, m~. ~ '7i?(l -|- A^). Thus MIA can work well for nearly degenerate squark 
masses and, in general, its reliability can be checked only a posteriori. However, for its simplicity, it has been widely 
used as a model i ndep endent analysis to find the constraints on different off-diagonal parts of squark mass matrices 
from experiments Q|. Therefore, we use MIA to calculate Wilson coefficients in the paper. 




As it is shown that both Br and CP asymmetries depend significantly on how to calculate hadronic matrix elements 
of local operators^. Recently, two groups, Li et al. [H Elj and BBNS [H [H, have made significant progress 
in calculating hadronic matrix elements of local operators relevant to charmless two-body nonleptonic decays of B 
mesons in the PQCD framework. The key point to apply PQCD is to prove that the factorization, the separation 
of the short-distance dynamics and long-distance dynamics, can be performed for those hadronic matrix elements. 
It has been shown that in the heavy quark limit (i.e., nib — > oo) such a separation is indeed valid and hadronic 
matrix elements can be expanded in as such that the tree level (i.e., the order) is the same as that in the naive 
factorization and the as corrections can be systematically calculated |l7j. Comparing with the naive factorization, 
to include the ag correction decreases significantly the hadronic uncertainties. In particular, the matrix elements of 
the chromomagnetic-dipole operators Qg^ have large uncertainties in the naive factorization calculation which lead to 
the significant uncertainty of the time dependent CP asymmetry in SUSY models [l^. The uncertainties are greatly 
decreased in BBNS approach p^. In the paper we shall use BBNS's approach ( QCD factorization ) to calculate the 
hadronic matrix element of operators relevant to the decays B — > (pKs, rj' Ks up to the order. 

The experimental results, S^k is negative and Sr/'K is positive but smaller than 0.7, which implies that new CP 
violating physics affects B 4>Ks in a dramatic way but gives B rj' Ks a relatively small effect, which has been 
thought to be problematic [2^. To solve the problem, Khalil and Kou invoke to have operators with opposite chirality 
because the decay constant of rj' is sensitive to the chirality of the quarks and that of 4) is independent of the chirality 
of the quarks, operators with opposite chirality give contributions with opposite signs to -B — > rj' Ks but contributions 
with the same sign to i? — > dKc; [lol| . In MSSM the LL {LR) and RR {RL) mass insertions give contributions to 
operators with opposite chirality respectively. It is shown in Ref.^3| that (1) it is possible to explain the experimental 
results if having both the LL and RR (or LR and RL) mass insertions simultaneously, and (2) it is impossible to 
explain the experimental results if having only the LL and/or LR (or RR and/or RL) insertions. The first claim 
is certainly valid in general. However, the second claim is of the consequence of using the naive factorization to 
calculate hadronic matrix elements of operators. The current-current operators contribute to _B — > rj' Ks but not 
to i? — > 4>Ks. Their effects are not significant in both the naive factorization approach and the QCD factorization 
approach. However, the as corrections of hadronic matrix elements of QCD penguin operators and chromomagnetic- 
dipole operators have significant effects and consequently make the claim (2) in Ref.[i3| not valid. We show in the 
paper that one can explain the experimental results, S^k is negative and Srj'K is positive but smaller than 0.7, with 
a fiavor non-diagonal mass insertion of any chirality when as corrections of hadronic matrix elements of operators are 
included. We show that in the case of LL and RR mass insertions the Higgs mediated contributions to S^k alone can 
provide a significant deviation from the SM in some region of parameters and a possible explanation of S,yK in la 
experimental bounds in a very small region of parameters, satisfying all the relevant experimental constraints. When 
including all the SUSY contributions, there are regions of parameters larger than those in the case of including only 
the Higgs mediated contributions in which the theoretical predictions for S^k and Sr/'K are in agreement with the 
data in la experimental bounds. In the case of LR and RL insertions a satisfied explanation can also be obtained. We 




* The chargino contributions liave been studied in Ref. I 111 
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show that the branching ratio oi B ^ v'^^S, which in SM with the naive factorization is smaller than the measured 
value, can be in agreement with data due to SUSY contributions in quite a large regions of parameters. 

The paper is organised as follows. In section II we give the effective Hamiltonian responsible for the b ^ s transition 
in MSSM. We give the Wilson coefficients of operators including those induced by NHBs at LO in MIA with double 
insertions. In Section III we present the decay amplitudes and the CP asymmetries S^k and Sn'K- In particular, 
the hadronic matrix elements of NHB induced operators to the as order are calculated. The Section IV is devoted to 
numerical results. We draw conclusions and discussions in Section V. The loop functions and some coefficients of the 
matrix element of the effective Hamiltonian as well as implementation of 77-77' mixing are given in Appendices. 



II. EFFECTIVE HAMILTONIAN 

The effective Hamiltonian for h ^ s transition can be expressed as0, |^ 

^ p—u,c ^ 2— 3,. ..,16 

+ C7^ Qy^ + Csg Qsg + C'j^ + C'^g Q'sg ) + ^.c. (5) 

Here Qi are quark and gluon operators and are given by ^ 

Ql = {SaPl3)v-AiPl3ba)v-A, Q2 = {SaPa)v-AiPl3bfj)v-A, 

Q3{5) = {Saba)v-AY^{qi3qi3)v-(+)A, Q4(6) = (sQ&/3)y-A ^(q/3<7a)y-(+)A, 

Q7(9) = ■:^iSaba)v-AY^eg{qpq0)Y_^_(^_)ji, (58(10) = ■7:i^ab(3)v~AY^eq{qi3qa)v+(-)A, 
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Qii(i3) = {sb)s+pY] —iqq)s-i+)p , 

q 

<3l2(14) = {Sibj)s+p'Y qi)s-( + )P ^ 

9 

Qi5 =sa^''(l + 75)6V ^ga^.(l + 75)g, 



Q16 = Si a^^il + 75) bj ^ —qj + 75) qi 

q "^^ 

g 



= ^m,s-„a'^''G;,^(l + 75)&/5, (6) 



where {qiq2)v±A = gi7^(l ± 75)92, iqiq2)s±p = ± 75)^2 p = u,c, g = u,d,s,c,b, is the electric charge 
number of q quark, Aq is the color SU(3) GcU-Mann matrix, a and (3 are color indices, and i^^^ (G^jy) are the photon 
(gluon) fields strength. 

The primed operators, the counterpart of the unprimed operators, are obtained by replacing the chirality in the 
corresponding unprimed operators with opposite ones. We calculate the SUSY contributions due to gluino box and 
penguin diagrams to the relevant Wilson coefficients at the niw scale in MIA with double insertions, as investigated 
in Ref. , which are non- negligible if the mixing between left-handed and right-handed sbottoms is large and results 
are 



t For the operators in SM we use the conventions in Ref. llSi where Qi and Q2 are exchanged each other with respect to the convention 
in most of papers. 

t Strictly speaking, the sum over q in expressions of Qi (i=ll,...,16) should be separated into two parts: one is for q=u, c, i.e., upper type 
quarks, the other for q=d, s, b, i.e., down type quarks, because the couplings of upper type quarks to NHBs are different from those 
of down type quarks. In the case of large tan/3 the former is suppressed by tan~^ /3 with respect to the latter and consequently can be 
neglected. Hereafter we use, e.g., C^^ to denote the Wilson coefficient of the operator Qn = {sb)<j^_p |^(cc)s_p. 
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where 



and X = TO? /to? with TO^i and ms being the common squark mass and ffluino mass 

respectively. The one-loop functions in Eq. Q are given in Appendix A. The parts of Wilson coefficients c['^ (i— 
3,. ..,6,77, 8g) which are obtained by single insertion are the same as those in Ref . |2^ . Differed from the single insertion 
results, the LR or RL insertion also generates the QCD penguin operators when one includes the double insertions. 

For the processes we are interested in this paper, the Wilson coefficients should run to the scale of 0{nrn,). Ci — Cio 
are expanded to 0{as) and NLO renormalization group equations (RGEs) should be used. However for the Cgg and 
C^j., LO results should be sufficient. The details of the running of these Wilson coefficients can be found in Ref. 
|2l] |. The one loop anomalous dimension matrices of the NHB induced operators can be divided into two distangled 
groups [2^ 1 



The operator of Qi \ is defined as Qi 
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Here and hereafter the factor, j^, is suppressed (i.e., the anomalous dimension matrix for (5ii,i2 is ^(f^^) ^ etc.). 
For Q[ operators we have 



(11) 



Because at present no NLO Wilson coefficients Cl , i=ll,...,16, are available we use the LO running of them in the 
paper. 

There is the mixing of the new operators induced by NHBs with the operators in SM. The leading order anomalous 
dimensions have been given in Refs.p5Ll2^. We list those relevant to our calculations in the following. Defining 
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The mixing of Qii,i2 onto the QCD penguin operators is 



12+1 




1,2, 








Oi 


-1/3 


1 


02 


-1 







28/3 


-4 




20/3 


-8 



(12) 



(13) 



JMQ) 





Qs 


Qi Qb 


Qe 


On 


1/9 


-1/3 1/9 


-1/3 


O12 












For the mixing among the primed operators, we have 



and 7 



(M'Q') 



7 



{MQ) 



(14) 



(15) 



The mixing of the new operators induced by NHBs with the operators in SM has non-negligible effects on the Wilson 
coefficients of the SM operators at the 0(mf,) scale. In particular, the Wilson coefhcient of the chromo-magnetic 
dipole operator Cgg at the 0{mi,) scale, which has a large effect to Smk {M — <j>, rj'), can significantly enhance due to 
the mixing. To see it explicitly we concentrate on the mixing of Oi (for its definition, see Ea. l|12|) ) onto Qsg- Solving 
RGEs, we have 



Oagi^^) = 

Mm) = 

V = 

where V and jaa are given by 



^(Aio)(r/(M)^"/"'° - 77(m)^^'^'/'^°) + Csgi^io)v""''^^''", 

c=l,...,4 

cc 7858(7/ 

a,6=l,...,4 



+ 2/3o/) V-^ = diag(7ii ,722,733 , 744]- 

with / being the 4x4 unit matrix. Using 
and Eq. (jTHJ), Eq. reduces to 

Csgi^^) = 0.68C8g(Mo) - 3.2Ci3(Aio), 

where Ci{fio) = ^Ci3(^o) has been used. 

In our numerical calculations we neglect the contributions of EW penguin operators Q7,...io since they are small 
compared with those of other operators. 



(16) 

(17) 
(18) 

(19) 

(20) 

(21) 
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III. THE DECAY AMPLITUDE AND CP ASYMMETRY 

We use the BBNS approach [TtL to calculate the hadronic matrix elements of operators. In the approach the 
hadronic matrix element of a operator in the heavy quark limit can be written as 

{MKs\Q\B) = (MKs\Q\B)f[l + J2rna':], (22) 

where {MKs\Q\B) f indicates the naive factorization result and M — (p, r{ . The second term in the square bracket 
indicates higher order corrections to the matrix elements . We calculate the hadronic matrix elements to the otg 
order in the paper. In order to see explicitly the effects of new operators in the MSSM we divide the decay amplitude 
into three parts. One has the same form as that in SM, the second is for primed counterparts of the SM operators, 
and the third is new which comes from the contributions of Higgs penguin induced operators. That is, we can write 
the decay amplitude for B MKs as 



(23) 
(24) 




C'i) for B^(j)Ks 
-CI) iorB^iKs 



A° and A" will be given in subsections. 

The time-dependent CP-asymmetry Smk is defined by 

aMK{t) = -Cmk cos{AMsot) + Smk sin(AMsoi), (25) 

where 



^ - i-\MK| 2ImAM_fs' 

^J^K = -, , |, f2 > Smk = -I , |, f2 • (.2oj 



Here \mk is defined as 



_ (q\ AjB^MKs) 

The ratio {q/p)B is nearly a pure phase. In SM Xmk — e'^'' +0(A^. As pointed out in Introduction, the MSSM can 
give a phase to the decay which we call 0^^^^. Then we have 

^ = ^.(2/3+0— )M ^5Mi^ = sin(2/3-f</.susY) (28) 
\A\ 

if the ratio -pj = 1. In general the ratio in the MSSM is not equal to one and consequently it has an effect on the 
value of Smk, as can be seen from Eq.(|2Sl. Thus the presence of the phases in the squark mass matrix can alter the 
value of Smk from the standard model prediction of Smk — ^™-'^Pj/4>k ~ 0-7. 

A. Bl^^Ks 

A° for B^ to the as order, in the heavy quark limit is given as|5ll27j 



A" = {cl,\sj^s\0){K\sY'b\B) X ^ VptV^ 



ps 

p—u,c 



as + + as - ^(ay + ag + a^p) 



(29) 



where a^'s have been given in Refs. |5l l27| and are listed in Appendix B. The hadronic matrix element of the vector 
current can be parameterized as (Klsjf'blB) = Ff^^{q'^){p'^g+p'^) + {F^^^{q^) - Ff^^ {q'^)){m% - m^-X/g^. 
For the matrix element of the vector current between the vacuum and (j), we have (^|s7^6|0) = m^f^e'l^. 
A""- for B'^ 4'Ks, to the as order, in the heavy quark limit is given asQ 



A^ = A'\Ci) + A'\C,~*C'i), 

TO, / 1 4m 



A-{C.) = {^\-s^^s\Q){K\-srh\B){-VtbV:s 



^neu ^ --ai2 + ai5 

TOfc y 2 TOf, 



(30) 
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fli's in Eq. (|30|l are 
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where Sc = ml /ml and 



Gl= —<!>4,{x), GF{s,x)= dt\n[s-xtt] 
Jo ^ Jo 

GF^is) = f\x^^GF{i 
Jo ^ 



[s — ie,x) 



(31) 



(32) 



(33) 



with X — 1 — X and $0(2;) — 6xx in the asymptotic limit of the leading-twist distribution amplitude. In Eq. (|31|l the 
expressions of V' and H^^^ can be found in the Appendix B. In calculations we have set m„ ,i = and neglected the 
terms which are proportional to ml/ml in E q.ll32ll . We have included only the leading twist contributions in Ea. H31|l . 
Here we have used the BBNS approach [TtLIISI to calculate the hadronic matrix elements of operators. 



B. ^ n'Ks 

A" for B ifKs, to the as order, in the heavy quark limit has been given in Ref.i28| as 
A" ^ VpbV;sm%A; 



p—u.c 



a; = F^-^' i 



SupaliKr^') + 2(a3(X?7') - a,{Kr^')) - ^i^, {al{Kr^') - al{Kri')) + -{-a,{Kr,') + a,{KT^')) 



alUK) + MA- {aliv'K) aUr,' K)) - -a%{v' K) 



(34) 



where F^^^(M=ri', K) is a i? ^ M form factor calculated at — m\ or m^^ w 0. The a^'s in Eq. H34|l have been 
given in Ref.[28j and are listed in Appendix B. 

We calculate A" for B rj'Ks, to the as order, in the heavy quark limit and the result is 
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(35) 



(a) 




(a') 





(b) 



(c) 





(b') 



(d) 




(d') 



FIG. 1: vertex corrections 



where decay constants J^'^,, are defined in Appendix C, = ™6 Rk/^, M^,' = i?^,, and a"'s and a^'s are defined 
as 
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(37) 



and y , HM1M2 could be found in Appendix B. In Eas. H3t)l37|l Vs,2 and Vt,2 come from vertex contributions, P's from 
penguin diagrams, and Hs^3 from liard-scattering contributions, wliicli will be shown in the following. In numerical 
calculations we set rriu.d = so that the terms which are proportional niq (q=u, d) in Eq. H35|) are neglected. 

• Vertex Contributions 

The vertex contributions come from Fig.lQJ. In the calculations of vertex corrections we need to distinguish 
two diagrams with different topologies showed in Fig. when we insert different operators. In the calculation 
of twist-3 contributions we find the infrared divergence cancels only if we include all four diagrams for each 
topology and use the asymptotic form of twist-3 distribution amplitude, i.e., 4>p{x) = l,4>cr{x) = 6x{l — x), 

which can be combined into — j- f p ^^.^^ [T^ l29l| . By a straightforward calculation, we obtain 
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FIG. 2: The two diflerent penguin contractions 



• Penguin Contributions 



In the calculations of penguin contributions we also need to consider the difference between two diagrams with 
different topologies showed in Fig.|(21). We obtain that the twist-2 parts of penguin contributions P^j^ 2: P^i^2 
are 
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Jo X 



(40) 



(41) 



Compared with twist-2 contributions, the twist-3 projection yields an additional factor of x which has appeared 
in penguin contributions proportional to G^^ and G|g . We therefore find 



p; 



M2,3 



^(hn^-GiO)) + (h.^-Gil) 
rub \3 fi / V3 M 



G13 



-8G- 



16 



_21n^-GF(l) 
M 

_21n!^-GF(l) 



4Gi 

-8G 



1^2,.^] -GFW 



16 



nib 



_21n^-GP(sc) 
Ai 



(42) 



r,Tl,EW 
M2,3 



-2 ^ X 



!^(lln!^-G(0)) + (iln!^-G(l) 
nib \S fi / \3 /i 



+ [Gi3 + AfcGi4] 



_21n^-GF(l) 
A^ 

mb 



8 [7V,Gi5 + G 



16j 



-l-21n!;;^)-GP(l) 



-4 [Gi5 + A^cGie] ( -2 In - GF(1) ) - 4 [G1, + iY.Gfg] 



nib 



_21n^-GP(se 
A^ 



43) 



with 



GF{s,) = / da: GF(s„x)0Af.(x). 
Jo 



(44) 
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FIG. 3; hard-scattering contributions 




o, 

FIG. 4: Two-gluon emission from a quark loop. A second diagram with the two gluons crossed is 
implied. 



• Hard-scattering Contributions 

Hard-scattering contributions come from Fig.@ which read 



S,3 



-UflMi (PB[0 a R ^ 3vmB (Pb{0 (PMi W (PKhpKV) 



-3(?i + v) ^J.Mi(l3BiO(t^Mip{u)(t)M2p{v) - {u-v) ^iMi <Pb{S) ^"^'1^ ^ (j^Ahpiv) 



+U^Mi (t>B{Cj4>Ahp{u) 



6 



(45) 



where 4>' means the derivative of (p. There is end-point singularity in Ea. (|45() which can be treated in the way 
given in Appendix B. 

Because rj' is a flavor singlet there are three extra contributions related to the gluon content of rj' . They come 
from: the b sqa amplitude, spectator scattering involving two gluons, and singlet weak annihilation. As 
analysed in Ref.|23], the singlet weak annihilation is suppressed by at least one power of h^Qco/^b in the heavy 
quark limit. So we need to calculate the contributions from the h sgg decay and spectator scattering which 
have not included in the results given above. In SM the relevant calculations have been carried out and results 
have been given in Rcf. 28] . We calculate them with insertions of new operators. In order to make the paper 
self-contained we also list the SM results given in Ref. p^. 

Two gluon mechanism 



• The b sgg amplitude 



The amplitude for b sgg comes from Fig. Q. We use the formula (16) in Ref. directly. For an operator 
O = {sTib){qT2q)-, where the T's denote arbitrary spinor and color matrices, one has (assuming the two gluons 
to be in a color-singlet configuration) 



where 



A{b^ sgg)\p,^^=-{usT,Ub){giqi)g{q2)\tT iT2A)\0) 



1 12mg "^'^ 



48to2 



+ 0(l/m^) 



(46) 



(47) 
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FIG. 5: Spectator-scattering contributions to the B — » Krj' decay amplitudes. The shaded blob 
represents the r]^'^ g*g* form factor. 



When T = V ± A, we obtain the contribution to the B KP decay amplitudes 
When r = 5 ± P, we obtain that the decay amplitudes of _B ~> KP are given as 

Ci4 



A charm 



If W)k-(l+75)6|i^(p)>{(Cn + f 



C 
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(48) 



(49) 



• Spectator Mechanism 



The Fig.© gives the spectator scattering contributions. By a straightforward calculation, we obtain the hard 
spectator-scattering contributions to the B — )■ Krj' decay amplitude using two-gluon mechanism: 



^spe 



Jo ? JO 2/ 



with 



^ln^ + ^-G(0,y-)-G(l,y) 



+ (^4 + Ce) 



1^ In - (n^ - 2) G(0, y) - G(s„ y) - G{l,y) 
3 /I 



-2G8f -I- Gi3 



^ f ^ In - G(0, y)) + f | ha ^ - G(l, y) 
nib \3 / V-J M 



_21n^-GF(l,y) 



-4GiE 



4-21n!;;^-GP(l,,) 



8Gi, 



-2G8^ff + Gi3 
-4Gi5 



rrib 



-1 - 21n — - GF{l,y) 



8G 



16 



-21n^-GF(l,y-) 



-l-21n— -GF(l,y) 



(51) 



(52) 



(53) 



with Sg — {mq/mi,)'^ , Uf — h. Obviously there is singularity arising from the end-point region in the Ea. H5U|) 
which makes us to consider the effect of fc^ in the end-point region. After involved the fc_L effect of spectator 
quark and turned to 6-space, the conjugate space of fc_L, Ea. H5U|) now reads 



/jspcc 



Cp fsfK ml r bdb f d^dy &>B (C, 
Jq Ja 



X \ - [KoibV^) - KoibV^)] [Piiy) {^K{y, b) + r^^piy, b)) + r^P^'^^^piy, b)] 
'bKiib^c) Koib^) - Koib^y 



+niB 



yrriB 



Pr^{y)^Kiy,b) 



(54) 
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where ^^Bi£,,b), ^K{y,b), ^p{y,b) are corrected distribution amplitude including Sudakov factor [Tsl [16^ and 
Kq^^i-) are modified Bessel function of order 0(1). Numerically we find that the contributions from the two gluon 
mechanism are negligible. 



IV. NUMERICAL RESULTS 



A. Parameters input 



In our numerical calculations the following values are needed: 



• Parameters related mixing of 77 — 77' 

We use the following input : 



msmb 

• Wolfenstein parameters 

We use the Wolfenstein parameters fitted by Ciuchini et al|3^: 

A = 0.819 ±0.040 A = 0.2237 ±0.0033, 
p = p(l - aV2) = 0.224 ± 0.038, p = 0.230 ± 0.039, 

fj = r]{l~ AV2) = 0.317 ± 0.040, = 0.325 ± 0.039, 

7 = (54.8 ±6.2)°, Vp^+V^ = 0.398 ±0.040. 



(56) 



= fg ml = 0.0025 GeV^ , = fs {2ra\ - ml) = 0.086 GeV^ , (55) 

= 39.3° ± 1.0° . 

• Lifetime, mass and decay constants 

r(50) = 1.56 X 10-1^5^ Mb = 5.28GeV, m^ = 4.2GeV, 

mc = 1.3GeV, = lOOMeV, /s = 0.190GeV, fx = 0.158GeV. 

• Chiral enhancement factors 

For the chiral enhancement factors for the pseudoscalar mesons, /ip = Hi^Le.^ -^^e take 

Rx±,o — R^± — 1.2, 
which are consistent with the values used in [T^ Isij , and 



(57) 



• Form factors 

In the paper we need two form factors: F^^^'(O) = 0.34 and F^^^ . Compared with these rather well studied 
form factors, the form factors for B rj'^'^ are poorly known, which has hindered theoretical predictions for B 
decays involving rj^'^ very much for a long time. We adopt the following parameterization for the form factors 
{P = T] 01 r]')m\: 

i.B^P(0).Fi£+F,^^^^±^, (58) 



where Fi and F2 both scale like (A/rrif,)'^/^ in the heavy-quark limit. In our numerical analysis we set Fi 




Fo^^'^(O) = 0.28 and take F2 = 



B. Constraints from experiments 

We impose two important constraints from B — > Xsj and Bg — > fj,~^fi~ . Considering the theoretical uncertainties, we 
take 2.0 x 10^"* < Br(i? Xgj) < 4.5 x 10""*, as generally analysed in literatures. Phenomenologically, Br(_B Xgj) 
directly constrains |C7-y(mb)p + |Cy^(mb)p at the leading order. Due to the strong enhancement factor mg/mf, 

associated with single 62^^^^^ insertion term in Cy'^(mb), Jjs^^^^' (^ 10^^) arc more severely constrained than 
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°Br(B^ (|)K)xlO^ ' Br(B^ (|)K)xlO^ ' Br(B^ (t)K)xlO^ 



FIG. 6: The correlation between ^^Ks and Br(_B 4>Ks) for the insertion of only one kind of chirality. (a) is for the LR 
insertion, (b) is for the LL insertion with only SM and NHB contributions included, and (c) is for the LL insertion with the all 
contributions included. Current Icr bounds are shown by the dashed lines. 



623'^^^^ ■ However, if the left-right mixing of scalar bottom quark 1533'^ is large (~ 0.5), 82^'"^^^ is constrained to 
be order of 10~^ since the double insertion term 5^^'^^^'^ 
Bs tJ-^ in MSSM is given as 



be order of 10 ^ since the double insertion term S^^'^^^'^ 5^^^^^*^ is also enhanced by mg/mi,. The branching ratio 



\CQ,imb) - C'^Mb) + 2m(C7io(mb) - CloCm,))^] (59) 

where fh = m^/niB^- In the middle and large tan/3 case the term proportional to (C m — C'(q) in Eq. H59|l can 
be neglected. The current experimental upper bound of Br(i3s /i+zi") is 2.6 x 10~^ To translate into the 

constraint on Cq^^^^, we have 



^\CQAmw) - C'q^S^w)V + \CQ,A^nw) - C'^Jniw^ < d (60) 

Because the bound constrains |Cq. — Cg. | (i=l, 2), ^ we can have values of |CqJ and \Cq. \ larger than those in 
constrained MSSM (CMSSM) with universal boundary conditions at the high scale and scenarios of the extended 
minimal flavor violation in MSSM 7] in which \Cq.\ is much smaller than |Cq. |. At the same time we require that 
predicted Br of _B ^ Xs/i+/i^ falls within 1 a experimental bounds. 

We also impose the current experimental lower bound AM^ > 14:Aps~^ |3^ and experimental upper bound Br(i? 
Xsq) < 9% [33. Because 5^3^'^^'' is constrained to be order of 10~^ by Br(B Xgj), their contribution to AMs 
is small. The dominant contribution to AM^ comes from insertion with both constructive and destructive 

effects compared with the SM contribution, where the too large destructive effect is ruled out, because SM prediction 
is only slightly above the present experiment lower bound. 

As pointed out in section II, due to the gluino-sbottom loop diagram contribution and the mixing of NHB induced 
operators onto the chromomagnetic dipole operator, the Wilson coefficients Cg^ can be large, which might lead to a 
too large Br of B ^ Xsg. So we need to impose the constraint from experimental upper bound Br(i3 Xgg) < 9%. 
A numerical analysis for C'g^g=0 has been performed in Ref.^]. We carry out a similar analysis by setting both Csg 
and Cgg non-zero. 



C. Numerical results for B — * <j)Ks 



In numerical analysis we fix nig — rriq — 400GeV and tan/3 = 30. We vary the NHB masses in the ranges of 
91GeV <mh< 135GeV, 91GeV < mn < 200GeV with m,, < mn and 200GeV < < 250GeV for the fixed mixing 
angle a = 0.6, 7r/2 of the CP even NHBs and scan 6^^^ in the range < 0.05 for A=B and 0.01 for A ^ B (A 

= L, R). 

Numerical results for B (pKs are shown in Figs. |H1E| where the correlation between S'^Xs and Br(i? (t>Ks) is 
plotted. Fig. IHlis for the insertion of only one kind of chirality. El[a)iEl!b) andlHI^c) correspond to the LR insertion, 
the LL insertion with only NHB and SM contributions included, and the LL insertion with the all contributions 
included, respectively. We find that for the LR insertion (Fig. EJa)) and the LL insertion (Fig. EJc)) there are 
regions of parameters where S^Ks f^tHs in Icr experimental bounds and Br is smaller than 1.6 x 10^^. In the case of 
the LL insertion with only NHB and SM contributions included (Fig. Elb)) S'^/f^ > 0.3 because C[ = 0, i=ll,...,16 



^ Cq^^ ^ are the Wilson coefficients of the operators Qi \ which are Higgs penguin induced in leptonic and semileptonic B decays and their 

definition can be found in Ref. 34j. By substituting the quark-Higgs vertex for the lepton-Higgs vertex it is straightforward to obtain 
Wilson coefficients relevant to hadronic B decays. 



14 



c 
CO 



-e- 



Br(B^ (t)K)xlO 



Br(B^ (t)K)xlO 



Br(B^ (t)K)xlO 



FIG. 7: The correlation between ^^Ks and Br(_B — * 4>Ks)- (a) is for the LR and RL insertions, (b) is for the LL and RR 
insertions with only SM and NHB contributions included, and (c) is for the LL and RR insertions with the all contributions 
included. Current la bounds are shown by the dashed lines. 




Br(B^' ri'K)xlO^ ° °Br(B^ r|'K)xlO-'* ° °Br(B^' ti'K)x10^ 



FIG. 8: The correlation between Sn'Ks ^^'^ Br(_B r)' Kg) for the insertion of only one kind of chirality. (a) is for the LR 
insertion, (b) is for the LL insertion with only SM and NHB contributions included, and (c) is for the LL insertion with the all 
contributions included. Current Icr bounds are shown by the dashed lines. 



and consequently Ci, i=ll,...,16 can not be as large as in the case with both the LL and RR insertions due to the 
constraint from Bg — > The similar result is obtained for only a RL insertion or a RR insertion. In Fig. 

[Zta)j[Ztb), and[7Jc) correspond to the LR and RL insertions, the LL and RR insertions with only NHB and SM 
contributions included, and the LL and RR insertions with the all contributions included, respectively. We find that 
there are regions of parameters where Sij,Ks f^-Us in la experimental bounds and Br is smaller than 1.6 x 10~^ in all 
three cases. The difference among the three cases is that the size of the regions is different: the smallest ones are for 
the LL and RR insertions with only NHB and SM contributions included, and the biggest ones for the LL and RR 
insertions with the all contributions included, as expected. In the cases of the LL and RR insertions and the LR and 
RL insertions Sij,Ks can reach —0.9 or so. Comparing |HIa) IHIc) with[7fa) [7fc), one can see that the case with both 
the LR and RL insertions (both the LL and RR insertions) has parameter regions with negative Scf,Ks larger than 
those in the case with the LR (LL) insertion. 



D. Numerical results for B — > r]' Ks 



Input parameters used in this subsection are the same as those in last subsection. Numerical results for B rj'Ks 
are shown in Figs. where the correlation between Sn'Ks and Br(i? rj'Ks) is plotted. The Fig. |Hlis for the 
insertion of only one kind of chirality. |SIa),|SJb), and[Hfc) correspond to the LR insertion, the LL insertion with only 
NHB and SM contributions included, and the LL insertion with the all contributions included, respectively. In all 
three cases Srj'Ks ^ and in the case of the LL insertion with only NHB and SM contributions included the minimal 
value of Srj'Ks is 0.5. The similar result is obtained for only a RL insertion or a RR insertion. In Fig. El Ela), |^b), 
andEIc) correspond to the LR and RL insertions, the LL and RR insertions with only NHB and SM contributions 
included, and the LL and RR insertions with the all contributions included, respectively. We find that there are 
regions of parameters where Srj'Ks i^^^s in la experimental bounds and Br is larger than 2.7 x 10~^ and smaller than 
11 X 10~^ in the first and third cases. For the first case, i.e., the case of the LR and RL insertions, there is most of 




Br(B^' ri'K)xlO^ ° °Br(B^ Ti'K)xlO-^ " °Br(B^ Ti'K)xlO^ 



FIG. 9: The correlation between Sn'Ks s^nd Br(i3 — > r)' Ks)- (a) is for the LR and RL insertions, (b) is for the LL and RR 
insertions with only SM and NHB contributions included, and (c) is for the LL and RR insertions with the all contributions 
included. Current la bounds are shown by the dashed lines. 
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FIG. 10: The correlation between S^irs ^^'^ ^<t>Ks for the insertion of only one kind of chirality. (a) is for the LR insertion, (b) 
is for the LL insertion with only SM and NHB contributions included, and (c) is for the LL insertion with the all contributions 
included. Current la bounds are shown by the dashed lines. 




^^fL "^m. ^^vi 

FIG. 11: The correlation between Sn'Ks ^^'^ S^ks- {^) is for the LR and RL insertions, (b) is for the LL and RR insertions 

with only SM and NHB contributions included, and (c) is for the LL and RR insertions with the all contributions included. 
Current la bounds are shown by the dashed lines. 



region of parameters where Sri'Ks is positive. In the third case, i.e., the case of the LL and RR insertions with the ah 
contributions included, the region which corresponds to positive Sri'Ks increases. For the second case, i.e., the case of 
the LL and RR insertions with only NHB and SM contributions included, Sj^'k > 0.5 in all regions of parameters, in 
contrast to the cases of the LR and RL insertions and the LL and RR insertions with the all contributions included. 
The reason is that the new contributions come mainly from the combination Ci3{mw) — C[^{mw) which is constrained 
by Bs IJ-^l-t - It is the same reason that make Fig. ^h) similar to Fig. |HIb). We would like to note that our 
results on Br(i? rj' Ks) can agree with the experimental measurement j35j 

Br{B Tj'K) = {5.8tli) x 10"^ (61) 

Indeed, it has been shown that the data can be explained within theoretical uncertainties in SM without introducing 
any new mechanism |28j|. 

In order to show explicitly there is the region of parameters in which Sr/'Ks is positive and S^Ks is negative for a 
set of values of parameters we plot the correlation between S'^/f^ and Sri'Ks ii^ Figs. ^|^] The Fig. ^|is devoted 
to the case with an insertion of only one kind of chirality. Fig. llUf a) is for the LR insertion. Fig. Iiur b) and Fig. 
llOf c) are for only NHB and SM contributions included and the all contributions included, respectively, with the LL 
insertion. One can see from Fig. Iiur a) that there is the region in which S,yks is positive and Si^Ks is negative and 
their minimal values are and —0.6, respectively. Fig. llOf b') shows that with only the LL insertion there is no region 
in which S.q>Ks is positive and S^Ks is negative if one includes only NHB and SM contributions. The reason is that 
C'li isim-w) = without a RR insertion (see, eg. (TJ) so that the Bi{Bs ) upper bound limits the size of 

C'li.ia, which leads to S^Ks ^ 0-3, as shown in Fig. El^b). When switching on all SUSY contributions among which 
the contribution coming from the chromo-magnetic dipole operator is dominant because its Wilson coefficient can be 
large due to including the double insertions (see, Eq. Q ), there is the region in which Srj'Ks is positive and S^Ks 
is negative and their minimal values are and —0.7, respectively, as shown in Fig. Iiurc). We would like to point 
out that our result in the case with only an LL insertion is different from that in Ref. jlOl| because we include the 
as corrections of hadronic matrix elements and the effects of NHB induced operators. We have also calculated the 
correlation between S,j,Ks and Sri'Ks i^ the case with only an RL insertion and an RR insertion respectively. In the 
case with only an RL insertion, although there is the region in which Sri'Ks is positive and Sfj^Xs is negative but there 
is no region in which both Sn'Ks and S^Ks are in agreement with data in the la bound. In the case with only an RR 
insertion, there are some points in the parameter space for which both S^'Hs and S^Ks are in agreement with data 
in the la bound. 

Fig. ^] is for (a) the LR and RL insertions, (b) only NHB and SM contributions included and (c) the all con- 
tributions included with the LL and RR insertions. We can see from the figure that the region exists for all three 
cases. For the case of the LR and RL insertions, there is a small region in which both Sri'Ks S^Ks are negative. 
For the case of the LL and RR insertions with only NHB and SM contributions included, the region of parameters in 
which Sri'K is positive and Scj,Ks is negative is smaller than that in the case of the LR and RL insertions and there 
is no region in which both S^'Ks and S^Ks are negative. There is almost no region in which Sr/'Ks is smaller than 
—0.5 for the case of the LL and RR insertions with the all contributions included and also for the case of LR and RL 
insertions. 
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The numerical results are obtained for rrig = mq=400 GeV. For smaller gluino and squark masses, the Wilson 

coefScient Cg^ becomes larger, which could make the effect on Smk and Br larger. However, indeed the effect is 

limited due to the constraint from B X^g. For fixed nig, the Wilson coefficient C^'g is not sensitive to the variation 
of the mass of squark in the range about from 100 GeV to 1.5 Tev. Therefore, the numerical results are not sensitive 
to the squark mass and would have a sizable change when the gluino mass decreases. For the very big gluino and 
squark masses (say, several TeV), SUSY effects drop, i.e., one reaches the decoupling limit. 

V. CONCLUSIONS AND DISCUSSIONS 

In summary we have calculated the Wilson coefficients at LO for the new operators which are induced by NHB 
penguins using the MIA with double insertions in the MSSM. We have calculated the as order hadronic matrix 
elements of the new operators for B 4>Ks and B rj'Ks in the QCD factorization approach. Using the Wilson 
coefficients and hadronic matrix elements obtained, we have calculated the time-dependent CP asymmetries S^k and 
SjjiK and branching ratios for the decays B — > (pKg and B — > rj'Ks- It is shown that in the reasonable region of 
parameters where the constraints from B^ — Bs mixing , T{b —^ 57), T(b —^ sg), T{b — > s/i"'"/i~), B — > fJ.'^fJ-~ are 
satisfied, the branching ratio of the decay for B — > (fiKs can be smaller than 1.6 x 10~^, and S^k^ can be negative. In 
some regions of parameters Scj^Ks can be as low as —0.9. The branching ratio and the time dependent CP asymmetry 
of the decay for B — s- rj' Ks can agree with experiments within la deviation in quite a large region of parameters. In 
particular, our result in the case with only an LL insertion or an RR insertion is different from that in Ref.|lOl| because 
we include the as corrections of hadronic matrix elements of operators and the effects of NHB induced operators. 

It is necessary to make a theoretical prediction in SM as precision as we can in order to give a firm ground for finding 
new physics. For the purpose, we calculate the twist-3 and weak annihilation contributions in SM using the method 
in Ref. [s^ by which there is no any phenomenological parameter introduced. The numerical results show that the 
annihilation contributions to Br are negligible, the twist-3 contributions to Br are also very small, smaller than one 
percent, and both the annihilation and twist-3 contributions to the time-dependent CP asymmetry are negligible. 
The conclusion remains in MSSM and we have neglected the annihilation contributions in numerical calculations. 

In conclusion, we have shown that the recent experimental measurements on the time-dependent CP asymmetry 
in _B ^ (j)Ks and B rj'Ks, which can not be explained in SM, can be explained in MSSM if there are flavor 
non-diagonal squark mass matrix elements of second and third generations whose size satisfies all relevant constraints 
from known experiments [B — > Xsj,Bs — > fi^ji^,B — > XsH^ jj,~ , B — > Xgg, AMs, etc.). Therefore, if the present 
experimental results remain in the future, it will signal the significant breakdown of the standard model and that 
MSSM is a possible candidate of new physics. 
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Note added 

When completing the paper we received the reference In Ref. the new bound of Br{Bs is given 

as 

Br{Bs -> jJ.^fi~) < 5.8 X lO"'^ at 90% confidence level. 

The new bound would give a more stringent constraint on the contributions of NHBs. We have carried out preliminary 
calculations. The preliminary results show that it is still possible that NHB contributions to Smk {M = 0, 77') alone 
can make a sizable even significant deviation from the SM. In other words, the qualitative conclusion in the paper is 
still valid if using the new bound. 



17 



References 



[1] B. Aubert et al, BABAR Collaboration, Phys. Rev. Lett. 89 (2002) 201802; K. Abe et al, Belle Collaboration, 
arXiv:liep-ex/0308036 

[2] Aubert et al. (BABAR Collaboration), liep-ex/0207070 T. Augshev, talk given at ICHEP 2002 (Belle Collaboration), 

BELLE-CONF-0232; K. Abe et al., BELLE-CONF-0201 hep-ex/0207098 
[3] The Belle Collaboration, K. Abe et al, liep-ex/0308035 ;BELLE-CONF-0344); the talk given by T. Browder at LP2003, 

http://conferences.fnal.gov/lp2003/program/S5/browder_s05-ungarbled.pdj 
[4] M. B. Causse, hep-ph/0207070 G. Hiller, Phys. Rev. D66 (2002) 071502;A. Datta, Phys. Rev. D66(2002) 071702; M. 

Raidal, Phys. Rev. Lett. 89(2002) 231803; K. Agashe and CD. Carone, hep-ph/0304229 B. Dutta, C.S. Kim, S. Oh, Phys. 

Rev. Lett. 90 (2003) 011801; J.-P. Lee, K.Y. Lee, hep-ph/0209290 Y.-L. Wu, Y.-F. Zhou, hep-ph/0403252 
[5] C.-S. Huang and S.-H. Zhu, Phys. Rev. D68 (2003) 114020 arXiv:hep-ph/0307354|^ 

[6] M. Ciuchini, L. Silvestrini, Phys. Rev. Lett. 89(2002) 231802; L. Silvestrini, hep-ph/0210031";talk contributed at ICHEP02); 
S. Khalil, E. Kou, Phys. Rev. D67 (2003) 055009; R. Harnik, D.T. Larson, H. Murayama, A. Pierce, hep-ph/0212180| 

A. Kundu and T. Mitra, hep-ph/0302123 R. Arnowitt, B. Dutta and B. Hu, Phys. Rev. D68 (2003) 075008; J. Hisano, 
Y. Shimizu, hep-ph/0308255 C. Dariescu, M. A. Dariescu, N.G. Deshpande and D. K. Ghosh, hep-ph/0308305 Y. 
Wang, hep-ph/0309290 V. Barger, C.-W. Chiang, P. Langacker, H.-S. Lee, hep-ph/0310073 S. Mishima, A. L Sanda, 
prep-ph/0311068 N. G. Deshpande, D. K. Ghosh, hep-ph/ 0311332, B. Dutta, C. S. Kim, S. Oh, G. Zhu, hep-ph/0312389, 
C.H. Chen, C.Q. Geng, hep-ph/0403188 ' 

G. L. Kane et al., hep-ph/0212092 Phys. Rev. Lett. 90 (2003) 141803. 

J.-F. Cheng, C.-S. Huang and X.-H. Wu, Phys. Lett. B 585 (2004) 287 < arXiv:hep-ph/0306086| . 

B. Aubert, et al., BABAR Collaboration, Phys. Rev. Lett. 91 (2003) 161801. 
S. Khalil, E. Kou, Phys. Rev. Lett. 91 (2003) 241602. 

S. Baek, Phys. Rev. D67 (2003) 096004;D. Chakraverty, E. Gabrielli, K. Huitu, and S. Khalil, Phys. Rev. D68 (2003) 
095004. 

C. -S. Huang and X.-H. Wu, Nucl. Phys. B657(2003) 304. 
L. J. Hall, V. A. Kostelecky, S. Raby, Nucl. Phys. B267 (1986) 415. 

For a review, see: F. Gabbiani, E. GabrieUi, A. Masiero, L. Silvestrini, Nucl. Phys. B477 (1996) 321. 

H. -n. Li and H. Yu, Phys. Rev. Lett. 74 (1995) 4388; H.-n. Li and T. Yeh, Phys. Rev. D56 (1997) 1615. 
Y. Y. Keum, H.-n. Li and A. I. Sanda, Phys. Lett. B504 (2001) 6; Phys. Rev. D63 (2001) 054008. 
M. Beneke et al., Phys. Rev. Lett. 83(1999) 1914; Nucl. Phys. B591(2000) 313. 

M. Beneke et al., Nucl. Phys. B6O6(2001) 245. ' 

R. Harnik, D.T. Larson, H. Murayama, A. Pierce, | hep-ph/021218'ol 
Y. Nir, Nucl. Phys. Proc. Suppl. 117 (2003) 111. 

G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68(1996) 1125 |arXiv:hep-ph/9512380| . 
L. Everett, G. L. Kane, S. Rigolin, L. -T. Wang, and T. T. Wang, JHEP 0201 (2002) 022:^ 
S.W. Baek, J.H. Jang, P. Ko and J.H. Park, Nucl. Phys. B609(2001) 442 . 

J.A. Bagger, K.T. Matchev and R.J. Zhang, Phys. Lett. B412(1997) 77; M. Ciuchini et al., Nucl. Phys. B523(1998) 501; 

C. -S. Huang and Q.-S. Yan, hep-ph/9906493 A.J. Buras, M. Misiak and J. Urban, Nucl.Phys. B586 (2000) 397. 

F. Borzumati, C. Greub, T. Hurth and D. Wyler, Phys. Rev. D62(2000) 075005 . 

G. Hiller, F. Kriiger, hep-ph/0310219 

X.-G. He, J.P. Ma and C.-Y. Wu, Phys. Rev. D63 (2001) 094004 ; H.-Y. Cheng and K.C. Yang, Phys. Rev. D64(2001) 
074004. 

M. Beneke and M. Neubert, Nucl. Phys. B651 (2003) 225. 
B.V. Geshkenbein and M.V. Terentev, Phys. Lett. 117(1982) 243 . 

R. Kaiser and H. Leutwyler, Proceedings of Workshop on Nonperturbative Methods in Quantum Field Theory, Adelaide, 
1998 hep-ph/9806336 ; 

R. Kaiser and H. Leutwyler, Eur. Phys. J. C 17 (2000) 623 hep-ph/0007101 . 
Hai-Yang Cheng, Yong-Yeon Keum and Kwei-Chou Yang, Phys.Rev.D65 (2002) 094023. 
M. Ciuchini et al, J. High Energy Phys. 07(2001) 13. 
CDF Collaboration, F. Abe et al., Phys. Rev. D57(1998) 3811. 

Y.-B. Dai, C.-S. Huang and H.-W. Huang, Phys. Lett. B390(1997) 257; C.-S. Huang and Q.-S. Yan, Phys. Lett. B442(1998) 
209 ; C.-S. Huang, W. Liao, and Q.-S. Yan, Phys. Rev. D59(1999) 011701. 
K. Hagiwara et al., Phys. Rev. D66 (2002) 010001. 

A. Stocchi, Nucl. Phys. Proc. Suppl. 117(2003) 145 arXiv:hep-ph/0211245 . 

A. Kagan, hep-ph/9806266 T.E. Coan et al. (CLEO Collaboration), Phys. Rev. Lett. 80 (1998) 1150. 
J.-F. Cheng and C.-S. Huang, Phys. Lett. B554(2003) 155 . In the paper the constraint from B is not imposed. 

R. Harnik, D. T. Larson, H. Murayama, A. Pierce, arXiv:hep-ph/0212180j 

D. Acosta et al. (CDF Collaboration), hep-ex/0403032 

T. Feldmann, P. KroU and B. Stech, Phys. Rev. D 58(1998) 114006 hep-ph/9802409|; 



T. Feldmann, Int. J. Mod. Phys. A 15(2000) 159 hep-ph/9907491 



18 



Appendix A Loop functions 

In this Appendix, we present the one-loop function of Wilson coefficients in this work. 



with 



^1,2(2; 
Pi, 2(2; 
^242,4,34(2; 
^1,2(2; 
^1,2(2; 
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dx 
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2 9x2 \y^- 
x" d^Ci,2{x,y) , 
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fi2{x) = 9fi{x) + f2{x), hi{x) = 9/3(2:) + /4(2:)/2 

where Bi 2(x,y) and Ci 2{x,y) are defined in Ref. 39] and /i, 2. 3,4.60 in Ref. [13]. 

Appendix B Coefficients ai 



(62) 



We shall give the explicit expressions of coefficients of the matrix element of the effective Hamiltonian in SM 
[T8LI28I ] which are not given in the content. For the integral which contains end-point singularity, we give a corrected 
one by including transverse momentum effects of partons and the Sudakov factor. The coefficients can generally 
be divided into two parts a,{MiM2) = ajj(MiM2) + aiji{MiM2): 
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where M2 means the "emission" meson, which is </) in the process B Kef) and is ?/ or K in the process B rj'K. 
The vertex contributions in Eq. H63|) are given as foUowing: 



pi 

12 hi— -18+ / dxg{x)^M2ix) 
M Jo 



12 In— -6 + / dxg{l - x)<S>M2{x) ; 



lix) = 3 



1 - 2x 
l-x 



In; 



2Li2(x) - In^a; 



2 hi a; 

l-x 



— (3 + 2i7r) lux — (x 1 — x) 



(63) 



The constants 18 and 6 are specific to the NDR scheme. Next, the penguin contributions are 



pp 



Ci 
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(64) 



where ny = 5 is the number of fight quark flavours, and s„ — 0, Sf. — (mc/mb)^ are mass ratios involved in the 
evaluation of the penguin diagrams. The smafi contributions from electroweak penguin operators are neglected in 
2 within our approximations. Similar comments apply to H67|l below. The function GAf2(s) is given by 



Gm2 (s) 



dx G(s — ie,l — x) (x) 



(65) 
(66) 



Compared with the twist-2 terms, the twist-3 terms do not have factor 1 — x, which is cancelled in the calculation. 
We therefore find 



G{s,x) = —4 / duu{l — u)\n[s ~ u(l ~ u)x] 
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(67) 



with 



Gm2 {s) = \ dx G{s -ie,l- x) $f ^ (x) . 
Jo 

The hard-scattering contributions HmiM2 ^i'g defined as following: 



HmiiP oc J d£,dudv 



We assume 



<^b{x) = Nbx^{1 - x)2exp 



2^1 



(68) 



(69) 



(70) 



with normalization factor Nb satisfying dx^^fa^) = I7 which is the popular used form in literature. Fitting various 
B decay data, lob is, determined to be 0.4 GeV 16]. With the Sudakov effects to cancel end-point singularity, Ea. H69|) 
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turns into 



r 
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Calculating the upper integral in b space, we obtain the hard spectator scattering contribution 
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(72) 



where Ki, li are modified Bessel functions of order i and ■ ^b , , are B, Mi, M2 corrected meson amplitudes 

with the exponentials 5 b , Si^[^ and Si^[^ respectively |0, ll^ 



Appendix C. Implementation of Tj—rj' mixing 

Though in the calculation of weak decay amplitudes we only consider the case with an 77' meson in the final state, 
we have to include the effect from rj-rj' mixing. In the following discussions we use the FKS scheme to describe the 
mixing between 77-77' [281 |4ll| . The matrix elements of local operators evaluated between the vacuum and 77^^'^ of the 
flavor-diagonal axial-vector and pseudoscalar current densities, read 



{Piq)\qri5q\0)= --^ .fp , 2m,{P{q)\qroq\0) = ^^P > 



(73) 



(P(g)|s7^75s|0)=-*/|.gr 



2ms{P{q)\sj5s\0) = -ih'p, 



where q — u or d. We assume exact isospin symmetry and identify vriq = \{mu + rrid). We also need the anomaly 
matrix elements 



{P{q)\^G^^,G^^n^)=ap, 



where we use the convention 
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-1) 



(74) 



(75) 



for the dual field-strength tensor. In all cases P = rj or rj' denotes the physical pseudoscalar meson state. 

The Ea. (|73|l and Ea. H74|) including ten non-perturbative parameters fp, hp, and ap, which however are not all 
independent, connect by 



d^.iqri5q) = 2tm, q^,q - ^ G% G"'^'' (76) 
(and similarly with q replaced by s) yields four relations between the various parameters, which can be summarised 



ap = 



p IfLp 



^/2 



= hp - fp nip . 



(77) 



which leaves us with six independent parameters now. 

The physical states are related to the flavor states in the FKS scheme by The 
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where \rjq) = {\uu) + \dd))/^/2 and |?7s) = |ss), then the same mixing angle appHes to the decay constants fp and /ip 
with the normaUzation given by (|73|) . We can re-define the decay constants in terms of fs,q and a mixing angle (j) as 

f^ = fqCOS(j), /^ = -/sSin0, 

(79 j 

and an analogous set of equations for the hp. Inserting these results into H77|l allows us to express all ten non- 
perturbative parameters in terms of the decay constants fq, fs and the mixing angle ((>. We obtain 

hq — fq (m^ cos^ 4> + TTi^' siu^ 0) — V2fs (m^, — m^) sine/) cos 4> , (80) 

hs — /s (m^, cos^ + TO^ sin^ (/i) ^ (m^/ — m^) sine/) cos (/>, (81) 

V2 



and 



^ V? '■■^'^"^^ ~ ^'^'^ ^ ''^ sin0cos0(-/5 sin(/) + \/2/s cos^) , (82) 

^ (/q m,j, — /ig) sin (/) = sin (j) cos </> (/^ cos + v 2/s sin (^) . (83) 

V 2 V 2 



